The Positronium state in quartz 



B. 



Barbiellini^ and P.M. Platzman-^ 
^ Department of Physics, Northeastern University, Boston, Massachusetts 02115, USA 
^ Bell Laboratories, Alcatel-Lucent, 600 Mountain Avenue, Murray Hill, New Jersey 07974-, USA 



o 
o 

Ph. 
< 



o 



> 
in 

On 
O 

o\ 
o 



X 



The positronium state in quartz is described by a linear superposition of two states: the first 
describing the free positron in the crystal and the second corresponding to a positronium Bloch 
wavefunction in the lattice. The condition for positronium formation in the electron gas is deduced 
by using variational calculations of the positron binding energy to the electron system. The self 
annihilation parameter k introduced in positron lifetime experiments can be properly justified by 
using the mixed state. A variational method to calculate « is proposed. 



PACS numbers: 36.10.Dr, 71.60.+Z, 78.70.Bj, 82.30.Gg 

Positronium (Ps), first observed in 1951 [l|, is a unique 
probe atom to test bound state QED Ps is formed 
in two spin states, 75% as ortho-Ps with parallel spins 
(triplet state) and 25% as para-Ps with antiparallel spins 
(singlet state). The energy difference between these spin 
states (hyperfine splitting) is about 8.4 x 10"'* eV. Para- 
Ps annihilates in vacuum mainly into photons of 511 keV. 
However, ortho-Ps annihilates in vacuum into three pho- 
tons in order to conserve the spin. The calculated life- 
time in vacuum for para-Ps is Tp = 124 ps while the 
corresponding value for ortho-Ps is Tq = 142 ns, which is 
in good agreement with recent experimental results 0]. 

Near surfaces the positron wave function overlaps with 
electrons outside the Ps. Therefore, the annihilation with 
electrons having an antiparallel spin decreases the ortho- 
Ps lifetime. This process is called pick-off annihilation 
and results in two photons. For instance, a quartz surface 
binds a ortho-Ps atom with an attractive van der Waals 
interaction [3] and the lifetime is then reduced to about 
10 ns H. 

Quartz is also one of the the rare materials, where Ps 
is delocalized in a Bloch state leading to peaks with neg- 
ligible widths at reciprocal lattice vectors in the angular 
correlation of annihilation radiation of the positron anni- 
hilation (AGAR) spectrum. Brandt et al [y] by studying 
the AGAR in quartz saw a spectacular fine structure. 
Greenberger et al. Q demonstated that this structure 
was the manifestation of Ps formation in the form of a 
Bloch state. With the help of high resolution two dimen- 
sional AGAR 8] it has been possible to clearly extract 
the square modulus of the Ps Bloch function at the T 
state in momentum space. 

Moreover, Saito and Hyodo @ showed the existence of 
two well defined spin states ortho- and para-Ps by mea- 
suring two different lifetime components and by intro- 
ducing a self annihilation parameter k. Inside the quartz 
crystal, the ortho-Ps lifetime becomes as low as 270 ps 
and the quantum mechanical interpretation of the Ps 
state inside the material becomes a controversial topic 
p^ . In fact, in some alkali metals the positron lifetime 
is above 400 ps (i.e. much larger than the ortho-Ps life- 
time in quartz) and the positron-electron pair correla- 
tion function is very similar to the case of Ps in vacuum. 



However, one does not observe any signatures of real Ps 
formation in alkali metals such as singlet-triplet lifetime 
splittings or Ps spikes in the AGAR spectra or 3-photon 
annihilation (Tlj . 

In the quartz crystal, the fine structure in momen- 
tum space also suggests a certain decoupling between 
the Ps and the electronic gas surrounding it. However, 
how could this condition be met inside the quartz crys- 
tal where the Ps sees a considerable amount of electron 
density? Moreover, how could the self annihilation pa- 
rameter K be different from zero without violating the 
quantum mechanical principle of identical particles? 

In this paper, we solve the problem of the quantum me- 
chanical interpretation of the Ps state in a solid by using 
a mixing of the state representing the positron freely 
moving in an electron system and the state ^'2 of a Ps at- 
tached to this system. This new state is characterized by 
a mixing angle 9, which can be calculated variationally. 
An interesting wavefunction ansatz ^"2 for the Ps bound 
state to metals has been given by Bergersen [12|]. Un- 
fortunately, the annihilation rates computed from this 
ansatz are much lower than the corresponding experi- 
mental values. However, we will see that \['2 may become 
stable when the electron density is less than the typical 
metallic values. 

We will start by studying the problem of a positron in 
a homogeneous electron gas. This system is completely 
characterized by the electron gas parameter 
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where n is the electron density. For instance, the Fermi 
momentum in atomic units is pp = 1.92/rs and the 
Thomas Fermi screening wave number is qtf = 1-56/y^ 
[1^. Variational Quantum Monte Garlo (QMG) calcula- 
tions can be performed by using simulation cells contain- 
ing N = 226 electrons and one positron [3, llBl- The 
variational wave function in the standard Jastrow-Slater 
form is given by 
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where is a Slater determinant formed by plane waves 
for spin up electrons, is a similar Slater determinant 
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for spin down electrons, ip+ is the positron wave function 
and J is a Jastrow factor used by Boronski The 
wave function 95+ is constant since the positron is in the 
lowest energy state. The Jastrow factor can be factorized 
as 



J — J eeJ epJ eepi 



(3) 



where Jee,Jep and Jeep contain the electron-electron, 
electron-positron and electron-electron-positron correla- 
tions respectively. The factor J has 12 variational pa- 
rameters, which have been optimized using the stochastic 
gradient approximation (SGA) [H, [l^. The positron- 
electron binding energy Eb shown in Fig. 1 is the nega- 
tive of the electron-positron correlation energy given by 



Ee =< ^lli/l^-l > - < *0|i^|«'0 >, 
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where H is the Hamiltonian of the system and 'I'o is the 
state of the system non interacting with the positron. 
Boronski by performing variational QMC calcula- 
tions for ranging from 2 to 10 has observed an anoma- 
lous behavior of the annihilation rates at Ts = 6. He 
suggested that this effect is an indication of Ps forma- 
tion. 

To check this hypothesis, we have considered the trial 
wavefunction proposed by Bergersen 
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where A is the antisymmetrization operator enforcing the 
Pauli principle, ^'at-i is the wavefunction of the system 
without an electron and (j){ri,rj) is the Ps wavefunction 
containing variational parameters. Similar wavefunctions 
have been used to describe excitons in solids [13 • Un- 
like the wavefunction used by Boronski, Bergersen's 
original ansatz neglects the role of the electron-electron 
correlation in the calculation of binding energy Eg. This 
energy can parametrized by 
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where E-p^ is the Ps binding enegy in vacuum and 77 = 
2.2. We have verified as shown in Fig. 1 that for > 
6 the positron binding enery Eb given by Bergersen is 
lar ger than Eb obtained from QMC calculations using 
^1 l3- Thus, this observation indicates that Ps becomes 
stable for > 6. Nevertheless, a more accurate study 
introducing a Jastrow factor in Bergersen's ansatz may 
change the critical value of r^. 

Next, we study the possibility in quartz to have a mix- 
ing 



1*3 >= cos(6')|^i > +sin(6l)|*2 > 
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The state ^'i is the same as *i except that the Kohn- 
Sham orbitals for quartz replace the planes waves of the 



homogeneous electron gas in the Slater determinants. 
Moreover, the corresponding positron state (p+ is given 
1.4 
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FIG. 1: Positron binding energy Eb in units of Eps for the 
Jastrow-Slater and Bergersen wavefunctions as a function of 
rs- The Jastrow-Slater points are obtained from the QMC 
calculation by Boronski 



by a periodic function with an energy / with respect to 
the electron conduction band. A value / = 1.7 eV is 
obtained using the positron and electron chemical poten- 
tials from density functional theory [l^ and the energy 
gap of Eg = 8.9 eV. Likewise, the state ^'2 is the same 
as 5*2 except that the corresponding ^'at-i is the quartz 
groundstate wavefunction and the corresponding is a 
Ps Bloch wavefunction. The variational parameters of 
*3 can be optimized by using the SGA [1^ [l^. The 
mixing angle and the self-annihilation parameter k @ 
are related via the formula 



sin2(6'). 



(8) 



The angle 6 or the parameter k can also be calculated 
with the SGA. This method can be used to verify the ex- 
perimental value K = 0.34 [9] and the Ps binding energy 
4.8 eV for quartz 

In conclusion, we have presented a state representa- 
tion of a Ps entangled with the electronic structure of 
quartz. This model is suitable to explain the fine struc- 
tures observed in the momentum density distribution and 
the singlet-triplet positron lifetime splitting. 
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